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A thin film of topological insulator (TI) on a dielectric substrate and a structure "bulk TI - dielec- 
tric film - bulk TI" are considered as natural double- well heterostructures suitable for a realization of 
counterflow superconductivity. The effect is connected with pairing of electrons and holes belonging 
to different surfaces of TI and transition of bound electron-hole pairs into superfluid state. The 
case of TI heterostructures subjected by a strong perpendicular magnetic field is considered. It is 
shown that such systems are characterized by two critical temperatures - a mean-field temperature 
of pairing and a much smaller temperature of superfluid transition. The dependence of critical 
temperatures on magnetic field is computed. The advantages of TI based structures in comparison 
with GaAs heterostructures as well as graphene based heterostructures are discussed. 

PACS numbers: 71.35.Ji; 73.20.-r; 73.22.Gk; 

I. INTRODUCTION 

During last two decades spontaneous interlayer phase coherence in quantum Hall bilayers was the subject of com- 
prehensive investigations. According to the theoretical predictions- - — , a double-layer electron system subjected by a 
strong perpendicular to the layers magnetic field at the total filling factor of Landau levels close to vt — 1 should 
demonstrate unusual transport behavior. Such a behavior is connected with an interlayer pairing of electrons and 
holes belonging to zero Landau level (formation of stable magnetoexcitons) and a transition of the magnetoexciton gas 
into a superfluid state. Superfluid state is expected to reveal itself in a flow of antiparallel electrical supercurrents in 
adjacent layers and vanishing of Hall voltage. Experimentally, the effect was investigated in AlGaAs heterostructures 
by a number of groups^—. A huge increase of the counterflow conductivity and strong lowering of Hall voltage was 
observed at temperature below 1 K. Nevertheless, in these experiments a state with infinite counterflow conductivity 
was not registered. It can be accounted for the presence of unbound vortices^—, but the question is still open. 

Typical magnetic field at which the effect is observed in AlGaAs heterostructures is about 2 T. At such field the 
condition for the filling factor vx = 1, as well as, the requirement for the magnetic length £ = ^Jhe/eB be larger or of 
order of the interlayer distance d is fulfilled. But this field does not provide smallness of the Coulomb energy E c = e 2 /el 
(e is the dielectric constant of the matrix) in comparison with the distance between Landau levels: huj c — ft 2 /m*£ 2 
(m* is the effective mass of carriers). In such a case the validity of the lowest Landau level approximation (commonly 
used in theoretical studies) is questionable. The opposite inequality E c < hio c can be achieved at larger magnetic 
field (smaller magnetic length), but smaller £ requires smaller interlayer distance. Under decrease of d the interlayer 
tunneling amplitude increases that is negative factor for the counterflow superconductivity-^—. 

Graphene systems open new perspectives for a realization of magnetoexciton superfluidity in bilayers^—. In 
graphene the relation between Coulomb energy and the distance between Landau levels does not depend on magnetic 
field. The zeroth Landau level in graphene is separated from the nearest positive and negative levels by the gap 
A£?oi = \/2hvF/l: where vf ~ 10 6 m/s is the Fermi velocity in graphene. The inequality E c < AEqi is equivalent 
to e > a e ff/y2, where a e ff — e 2 /hvp ps 2.2 is the effective fine structure constant for the suspended graphene. 
The latter is fulfilled in graphene-based heterostructures, where Si02, AI2O3 or BN compounds are used as dielectric 
partsS&Si. 

As was shown in Refs. Il9ll20l . an imbalance of filling factors of graphene layers is required for a realization of 
magnetoexciton superfluidity in bilayer graphene structures. It is connected with an additional 4-fold degeneracy of 
Landau levels due to spin and valley degrees of freedom. This behavior is in similarity with one for i>i — 2 quantum 
Hall bilayers^,. The imbalance can be created by an electrostatic field applied perpendicular to graphene layers. 
Electrical field should be tuned simultaneously with magnetic field. It is required to keep the ratio E/B close to the 
value E/B = ac/e, where a ps 1/137 is the fine structure constant. For instance, for e = 4 and B pa 1 T electrical 
field E ps 5 ■ 10 3 V/cm is required. 

Discovery of topological insulators (see Refs. (26jj27| and references therein) stimulates new proposals toward real- 
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ization of superfluidity of spatially indirect excitons^r— . The ideas 2 ^— are based on that the surface of TI works 
as a natural two-dimensional conductor, while the interior of TI as a dielectric. The electron spectrum of TI surface 
states is similar to graphene spectrum - it contains Dirac cones. Therefore, TI surface should demonstrate the same 
quantum Hall behavior as graphene. In the TI surface the number of Dirac cones is odd, in particular, the most 
studied three-dimensional TI Bi2Se3 belongs to a so-called one-cone family. Therefore, TI systems have common with 
graphene systems advantage (smallness of Coulomb energy in comparison with the distance between Landau levels), 
but, at the same time, free from the disadvantage, caused by the additional degeneracy of Landau levels in graphene. 

In this paper we analyze two types of structures: a TI film on a dielectric substrate, and a heterostructure "bulk 
TI - dielectric film - bulk TI". In Sec. [IT] we obtain the zero-temperature phase diagram in the coordinates "the ratio 
d/£ - dielectric constants". It is shown that for the TI film on a substrate structure the diapason of allowed d/l is 
restricted from below and from above and it shrinks under increase of the dielectric constant of the substrate e s . For 
known TI dielectric constant eti is large (~ 10 2 ) and typical situation corresponds to e s /sti ?S 0.1. In this case the 
state with electron-hole pairing is expected in the diapason d/l = 0.5 2 (where d is the TI film thickness). For the 
" TI - dielectric film - TI" structure the diapason of allowed d/£ (d is the dielectric film thickness) is restricted only 
from above and it shrinks under decrease of the dielectric constant of the dielectric film Ed- For Ed/en ^ 0.1 the 
upper boundary for d/£ is small (< 0.2). 

In Sec. IIIII the mean-field critical temperature T m t and the temperature of the superffuid transition T s versus 
magnetic field is computed. The temperature T m f is given by the self-consistence equation for the order parameter. 
The temperature T s is the Berezinskii-Kosterlitz-Thouless transition temperature. We find that the strong inequality 
T m f -C T s is fulfilled. In TI film structures superfluid state is reached at temperatures higher than in "TI - dielectric 
film - TI" structures. The temperature T s is non-monotonic function of magnetic field and increases under decrease 
of e s and under increase of Ed- 



II. ZERO-TEMPERATURE PHASE DIAGRAM 



Let us consider electron surface states of a one-cone TI in a quantizing magnetic field directed perpendicular to the 
surface. The low-energy Hamiltonian has the form Hq = ±VF(p x a y — p y cr x ) + ma z , where Oi are the Pauli matrixes 
that act in the spin space, p = —iKS7 + eA/c is the momentum operator, A is the vector potential, m is the Zeeman 
splitting, and vp is the Fermi velocity that is the material parameter (typically, vf ~ 5 • 10 5 m/s). The eigenproblcm 
is given by the Dirac equation 

T**(p- + §A_)\/* t \ /%\ 
±iv F (p+ + %A + ) -m ){ * 4 J * { ttj. J ' {L> 

where p± = —ih(d x ± id y ) and A± — A x ± iA y . The upper (lower) sign in Eq. (TTJ) corresponds to the top (bottom) 
surface. 

In the Landau gauge A = (0, xB, 0) the eigenproblem (TTJ has the following solution. 
1. Zero Landau level. 

E = -m, <l>o,x(r)= e ( : ) , (2) 



where kx = X/£ 2 is the guiding center of the orbit, and L y is the size of the system in y direction. 
2. Nonzero levels 



E N = sign(iV) 

Tr 1 / 4 ^2\ N \\N\lL y £ y H^Zmg 



p-ikxVp — 1 

= ' °„ r „ j= i . I ^n^T^xT' ) - (3) 




where N = ±1,±2, . . .. 

In moderate magnetic field the Zeeman splitting m = gfisB/2 (where [ib is the Bohr magneton, and g is the 
gyromagnetic ratio) is small in comparison with the distance between Landau levels (mjE\ sa 3 • 1CP 3 at B = 1 T). 
Therefore, one can put m = in Eqs. Up to the phase factors, it yields the same eigenfunctions as for graphene. 

We imply that at zero magnetic field the Fermi level is tuned to the Dirac point E = 0. It can be done by the 
appropriate doping of TI^. In magnetic field the zero Landau level is shifted from zero. But if the number of surface 
carriers is conserved, the Fermi level is shifted as well and it coincides with zero Landau level. 
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The Hamiltonian of Coulomb interaction of electrons in zero Landau level reads as 

H =\Y,j d^r'V^dv-T'D-.piir^ir'):, (4) 

where Vi^ (r) is the Coulomb interaction between electrons located on i and i' working surfaces, 

P*(r)= *o,^i( r )*o,x 2 (r)c+ XiCi ,x 2 (5) 

X\ ,X2 

is the electron density operator in the second quantization representation, cf x (c^x) is the creation (annihilation) 

operator for electron in zero Landau level in the surface i, and : O : means the normal ordering of operator O. 
Substituting ([2|), ([5]) into (j4} we obtain the following expression for the Coulomb interaction Hamiltonian 

ii' X,X',q * z 

where Vn> (q) are the Fourier-components of the Coulomb interaction and S is the area of the system. 

In what follows we neglect the influence of outer boundaries on the interaction between working surface electrons 
and consider a model heretostructure "infinitely thick dielectric 1 - working surface 1 - dielectric 2 of thickness d - 
working surface 2 - infinitely thick dielectric 3". In the general case dielectrics 1,2, and 3 are characterized by different 
dielectric constants ei, £2, and £3, correspondingly. For such a structure the quantities Va>{q) read as 

(7) 
(8) 
(9) 
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The pairing of electrons from the surface 1 with holes from the surface 2 is characterized by the order parameter 

A x = (*\cf x C2 X m, (10) 

where | Vf) is the many-particle wave function. In ([TU)) the relation between the electron annihilation and hole creation 
operator Cix = hf x is taken into account. We consider the many-particle wave function 

l*> = II ( uc tx h tx + «) Ivac) = J] (uc+ x + vc+ x ) |0) (11) 

X X 

that is an analog of the wave function introduced in the Bardeen-Cooper-Schricffcr (BCS) theory of superconductivity. 
The u — v coefficients satisfy the relation \u\ 2 + \v\ 2 = 1 and are parameterized as u = cos(#o/2), v — e lVa sin($o/2), 
|0) is the state with empty zero Landau level, and |vac) is a "vacuum" state defined as |vac) = Yix c 2x\®)- 
The energy of the state (dU E = (ty\H\^) reads as 



E S 



Wcos 2 e Q - ( Jn ~t J22 (1 + cos 2 6 ) + (J11 - J22) cos^o I - Ji 2 sin 2 d 



8ir£ 2 
where 



2ire 2 q ^o\ 2 " W J £ 2 ^ 2 

is the energy (per particle) of the direct Coulomb interaction, and 



(12) 



w = ^( VMq) : v > M -v^)) = ^ m 



2tt 

are the energies of the intralayer and interlayer exchange interaction 



Jik = 7> / QV ik (q)e-^-dq (14) 
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In the state (fTTj) the filling factors of zero Landau level in the surfaces 1 and 2 are 

1± cos 6» , . 

^i(2) = 2 ■ ( 15 > 

They satisfy the condition v\ + v% — 1, but their difference v = v\ — v 2 = cos6*o (the filling factor imbalance) is 
a variable that should be found from the condition of minimum of the energy (|12p . and it is the function of the 
parameters of the heterostructure. 



The minimum is reached at 



1 ' dl 2(W+J 12 )-J 11 -J 22 1 < 



cos^o = •(-!, _ at 2{w+ J jUX_j^ < -1; (16) 

2(w+ J &)-j^-Jn > otherwise. 

Direct computation shows that 2(W + J12) — J11 — J22 > for any d and £j. 

In the state (JT]) the modulus of the order parameter is equal to |A| = sin6*o/2. Zero imbalance (cos#o = 0) corre- 
sponds to maximum order parameter. This case is realized in symmetric heterostructures (ei = £3). In asymmetric 
heterostructures the imbalance is nonzero that results in lowering of the order parameter. If the imbalance becomes 
maximum (9 = ±1) the order parameter goes to zero. Direct evaluation of (|16[) shows that the imbalance increases 
under decrease of d/i and it reaches the maximum at certain nonzero value of that ratio. Thus, there is the critical 
d/i below which the electron- hole pairing does not occur in the asymmetric system. 

A restriction on the parameter dj I also comes from the dynamical stability condition (the condition for the collective 
mode spectrum be real valued). To obtain the spectrum of excitations we follow the approach Ref. [34| and consider the 
many-particle wave function that takes into account fluctuations of the phase and the modulus of the order parameter 

I*) =11 ( cos 9 -f c tx + z lVx 6 f4x) |0>. (17) 
In the quadratic approximation the energy of fluctuations can be presented in the diagonal form 



E fl = lj2( m *M ) K ^ ( ™{q) ) ' ( 18 ) 

q ^ ' 

where 



m z (q) = -J—Y,(co S x -cos6 o )e-*i x (19) 



2V S x 



and 



9ttP 2 

V(q) = J^ r J2& x -Vo)e- iqX (20) 



S x 



are the Fourier-components of fluctuations of the filling factor imbalance and the phase of the order parameter. The 
matrix K in (|18p is 



K^{q) 
with the components 

K zz (q) = 2 (H(q) - F s (q) + F D (0) + cot 2 9 E(q)) , (22) 



1 



where 



K 9V ,(q) = - S m J 6 Z(q), (23) 



1 fV 11 (q) + V 22 (q) 



H (l) = ^2 ( o ~ V M ) ( 24 ) 
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E(q) = F D (Q) - F D (q), (25) 
F s(<?) = ^ y pJb&*tf 2 ) (K u (p) + V 22 (p)) e-^dp, (26) 

and 

1 pOO 2 2 

= 5- / pJ (pqt 2 )V 12 { P )e-^dp. (27) 
27r Jo 

In ([2l) |) .(|2T] ) J (g) is the Bessel function. 

The quantities m z (q) and (p(q) are canonically conjugated variables. They satisfy the equations of motion h(fi(q) = 
K zz (q)m z (q), Tirn z (q) = —K ipip (q)ip(q) from which the following expression for the collective mode spectrum comes 
from 

n(q) = y/K zz (q)K w (q). (28) 

Analysis of (|28p shows that at d/t larger than some critical one the quantity £l(q) becomes imaginary valued at a 
finite q. It signals for a dynamical instability of the system. 

For the symmetric structure " TI - dielectric him - TP the range of existence of the superfluid state is given by the 
inequality d/t < d c \ 1 where d c \ is the function of the ratio ed/sri- It is shown in Fig. [T^,. On can see that the diapason 
of djt suitable for a realization of magnetoexciton superfluidity broaden under increase of Ed- For the asymmetric 
structure "TI him on a dielectric substrate" the range of existence of the superfluid state is d c2 < d/t < d c %, where 
d c i and d c2 are the functions of two ratios: £1/62 and £2/63. Here we specify £1 = 1, £2 = £ti = 80 and present 
dependencies of d c i and d c2 on the dielectric constant of the substrate £ s (Fig. [TJd). One can see that in some sense 
the situation is opposite to one for the "TI - dielectric film - TI" structure - the use of dielectric substrates with 
larger e shrinks the range of allowed d/t. We emphasize that filling factor imbalance can be controlled by a gate 
voltage applied to the system. Therefore, in systems subjected by an electrical field directed perpendicular to working 
surfaces the lower and upper critical d/t will differ from ones presented in Fig. [T] 



III. FINITE TEMPERATURE PROPERTIES 



Let us consider finite temperature behavior of the system in the framework of the mean-field approac h 35 i 36 . 
One obtains from (j6j) the following mean-field Hamiltonian 



H 



MF 



E 

X 



eicf x c iX - {Ji 2 Ac+ x c 2X + h.c.) 



where A = (cj^cix) is the mean-field order parameter, 
with vi = {cf x Cix) t the mean-field filling factors, 



f 



„ lim 

2-Kt 1 9->0 



Vui'l) I 'I ^) + V 12 (q) \V2-\ 



(29) 



(30) 



(31) 



Do = 



1 



„ lim 

2nt 2 g^o 



V 22 (q) 



v 2 

2 



V 12 (q) 



V\ 

2 



(32) 



and fi is the chemical potential. The quantities (|31l) . Q32p account the interaction with positive background as well. 
The Hamiltonian (|29[) is diagonalized by the u — v transformation 



dx = ua x 



■v*f3 + 



c 2 x = u*f3 x - vax, 



(33) 
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FIG. 1: Zero-temperature phase diagrams for the TI- dielectric film - TI heterostructure (a) and for TI film on a dielectric 
substrate heterostructure(b). 



where u — cos(0/2) and v = sin(6/2)e l¥ '. The condition of vanishing of non-diagonal terms in the transformed 
Hamiltonian yields the following relations 



. Ji 2 |A| 
sin 6= — , cos 6 



where 



ei - e 2 1 



W - 



•hi + J22 \ ~ J11 — J: 



22 



2 2 

The transformed Hamiltonian has the form 

Hmf = ^2 (E a a^a x + EpP^px) 



x 



with the spectrum 



E, 



(0) 



?2 + J2 2 |A|2± 



ei + £2 



The condition x/\ + V2 = 1 yields the relation 

1 = l + N F (E a )-N F (E fi ), 
where Np(E) is the Fermi distribution function. Relations (l37l) . (|38)) lead to the condition 

(-1 + (-2 = 0. 



(34) 

(35) 
(36) 

(37) 

(38) 
(39) 
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Note that Eq. ([39]) can be satisfied by the appropriate choice on the chemical potential /i. It corresponds to that the 
chemical potential is determined by the relation v\ + = 1. Under accounting Eq. Q39p the spectrum (|37p is reduced 
to £ Q = Ep = £ = v/? 2 + J? 2 |A| 2 . 

Eq. (|38|l is the first one in the set of three self-consistence equations. Two other equations read as 

e , £ 

A = ^tanhA. (40) 
2£ 2T y ' 

It follows from (|40p that at A ^ the following relation is fulfilled 

J 12 i> + 2£ = 0. (41) 

Eq. (gU) yields 

Jn — J22 



2(W + J12) - J11 - J22 



(42) 



that coincides with Eq. (|16[) under assumption that | Jn — J22I < 2(W + J12) — J11 — J22- If the latter inequality is 
not fulfilled, Eq (|4"Tj) cannot be fulfilled as well, and the order parameter A = (electron-hole pairing does not occur). 
Thus, in the paired state the filling factor imbalance does not depend on temperature and is given by Eq. (|42[) . 
Eqs. (|40[) yield the mean-field critical temperature of pairing^ 

T J12 |£| ,,„s 



2 ln^ 



PI 



where £ is given by Eq. (14"2"1) . 

The temperature T m / is the function of d/t. For a given sample the distance d is fixed and the parameter d/t 
depends only on magnetic field. Therefore, it is instructive to present the critical temperature as the function of 
magnetic field B. We choose Bd = fa/nd 2 units for B, where 0o = frc/2e is magnetic flux quantum (B/Bj = {d/t) 2 ). 
The quantity e 2 /d is used as the energy unit. In computation T m f we account the dynamical stability condition, 
implying that T m f = at d/t > d c \. The result of computations for two types of heterostructures is presented in 

Figi ' ^ 

The critical temperature T m f is not a temperature of superfluid transition. Superfluid transition temperature is 
given by the Kosterlitz-Thouless equation 

T s = ^ Ps (T s ), (44) 

where p s {T) is the superfluid stiffness. Equation (|44|) can be applied under assumption that the gas of bound electron- 
hole pairs exists. The latter requires T s < T m f. Evaluation of T s shows (see below) that, actually, the strong inequality 
T s <C T m f is fulfilled. Therefore, the temperature T m f can be interpreted as an analog of the ionization temperature. 

The superfluid stiffness is the coefficient of the expansion of free energy in the gradient of the phase of the superfluid 
order parameter: 

F = F a + ^ Ps (V^) 2 . 
We compute p a {T) as follows. We consider the many particle wave function 



X 



Eq. (1451) describes the state with nonzero space-independent gradient of the phase of the order parameter Vip = Q = 
(Qx, Qy)- To see that we neglect for a moment the fluctuations (Ox = 0q and ipx = 0) and define the space-dependent 
order parameter 

A(r) = rx{*)4>x>(v)mc+ x c 2X ,\^). (46) 

X,X' 
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FIG. 2: The dependence of the mean-field critical temperature on magnetic field for the TI- dielectric film - TI heterostructure 
(a) and for the TI film on a dielectric substrate heterostructure(b). 



In pj ip x (r) = {y/ir^Lyl)- 1 exp(-iXy/£ 2 - (x - X) 2 /2£ 2 ) is the one-particle wave function for zero Landau level. 
Direct calculation yields 



A(r) = 



sinflp iQ . r 



Eq. (147)) shows that Q is indeed the gradient of the phase of the order parameter. 
The free energy is given by the formula 



F = E {Q) -TS = E (Q) + T^ln (l - e- 1 ^ 22 ) 



where 



E (Q) = Eo(0) 



S 



sin 2 9 [F D (Q) - F D (Q)} 



(47) 



(48) 



(49) 



is the energy of the state (l4"5j) . S is the entropy of the gas of elementary collective excitations, f2(q, Q) is the spectrum 
of excitations, and Eo(0) is the energy given by Eq. (| 12|) . 

At Q = the spectrum i7(q, 0) = £l(q) (Eq. t\'28\i ) is isotropic. Anisotropy of f2(q, Q) i s connected with the 
appearance a specific direction in a system (the direction of the phase gradient). 

Considering the state (|45)l and repeating the same steps as in obtaining the spectrum (l28|) we finds 0(q, Q) for q 
directed parallel to the x axis. 

The energy of fluctuations has the form (|18p . where the matrix K depends on Q: 



K(q,Q) 



K zz (q,Q) K zv (q,Q) 



with q = qi x , 



K* ztp (q,Q) K w (q,Q) 
K zz (q, Q) = 2 (H(q, Q) - F s (q) + F D (Q) + cot 2 9 ~(q, Q)) , 



(50) 



(51) 
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K v Jq., Q) = i sin 2 #oS(q, Q), 



(52) 



^(q, Q) = iK zlp (q, Q) = i cosd Q [F D (\q_ - Q|) - F^flq + Q|)] /2, 



(53) 



ff(q, Q) = ^ ( Mgl+Mg) _ Fl2(?) cos( | q x Q|) ) e"^, 



and 



s (q> Q) = Fd{Q) - 



F^dq + QD+Fodq-Ql) 



The spectrum has the form 



fi(q, Q) = t/jr,,(q, Q)# w (q, Q) + tf, v (q, Q). 



(54) 



(55) 



(56) 



One can see that the expression (|56| is invariant with respect to rotation of coordinate axes. Therefore, the restriction 
q = qi x can be omitted and Eq. (|56l) yields the spectrum at the general q. The spectrum Eq. (|55|) is anisotropic 
since it depends on the angle between q and Q. 

Expanding (|48p in Q we arrive at the following expression for the superfluid stiffness 



Ps = psO + 5p s {T), 



where 



PsO 



sin 2 6 a 
8ir£ 2 



f d (Q)\q=o 



sin 2 9 Q 



P 3 V 12 (p)e ~dp 



is the zero-temperature superfluid stiffness, and 



, d 2 Sl(q,Q) 1 ^»(g,Q) 



(57) 



(58) 



(59) 



2=o 



is its temperature correction. In (|59p Nb(q) — (e^ q '' T — 1) _1 is the Bose distribution function. 

One can show that 8p s (T) < 0. We note that Eq. (l59l) generalizes the expression for the superfluid stiffness^. The 
case Ref. [53 corresponds to a Bose gas in a free space. In that case the Galilean transformation yields the following 
spectrum of excitations f2(q) = flo(s) + ' Q/M with ilo(q), the excitation spectrum in the reference frame where 
the Bose gas is at rest, and M is the mass of Bose particles. Then, Eq. (jBTJf is reduced to the common one 37 . The 
dependence of the spectrum f2(q, Q) Eq. (l56|) on Q is more complicate, and the first term in (|59| should be taken 
into account. 

The dependence of superfluid transition temperature on magnetic field is shown in Fig. [3] One can see that at all 
B the inequality T s -C T m f is satisfied. The temperature T s is a non-monotonic function of magnetic field. Superfluid 
state can be realized in the whole diapason of allowed d/£ {d c \ < d/i < d C 2), but at d/l near the upper or lower 
boundary of this diapason the temperature T s approaches to zero. The maximum temperature is reached in the 
middle of that diapason. 

For the TI film on a substrate structure higher temperature of superfluid transition correspond to lower dielectric 
constant of the substrate e s . For the " TI - dielectric film - TI" structure decrease of the dielectric film parameter 
results in lowering of the superfluid transition temperature. 



IV. CONCLUSION 



In conclusion, we have shown that TI heterostructures are suitable for realization of superfluidity of spatially indirect 
magnetoexcitons. The structure "TI film on a substrate" is preferable to the "TI - dielectric film - TI" structure. 
The main disadvantage of the latter ones is that for typical materials the thickness of dielectric layer that separates 
two TI should be rather small. For instance, for B = 1 T, = 12 and eti = 80 the film should not be thicker than 
8 nm. For the same parameters the thickness of TI film on a substrate with e s = 12 can be up to 50 nm. The TI 
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FIG. 3: The dependence of the superfluid transition temperature on magnetic field for the TI- dielectric film - TI heterostructure 
(a) and for the TI film on a dielectric substrate heterostructure(b). 



film structure has the problem of shorting of two working surface through the side surface. But this problem can 
be resolved by depositing of magnetic insulator on the side surface. It opens a gap in the side surface state energy 
spectrum: E = ± yj m 2 + Vpp 2 (with m = JrS z 3> T s , where Jh is the exchange energy with magnetic ions with spin 
S z ) and prevents interlayer leakage of the counterflow current through the side surface. 

Taking d = 10 nm we evaluate maximum temperature of superfluid transition T s about 1 K for TI film on a substrate 
structures and about 0.2 K for TI - dielectric film - TI structures. Graphene heterostructures are characterized by 
a little but higher— temperature of transition of magnetoexciton gas into superfluid state, but the disadvantage of 
graphene heterostructures is that magnetoexciton superfluidity can be realized only under application of interlayer 
gate voltage. 
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